Atmospheric Neutrino Fluxes 



oo 



Thomas K. Gaisser 
Bartol Research Institute, University of Delaware 
Newark, DE 19716 USA 



in 
o 
o 

<N ■ February 2, 2008 



Abstract 



Starting with an historical review, I summarize the status of calculations of 
| the flux of atmospheric neutrinos and how they compare to the measurements. 

OO 
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<N 1 Introduction 

o ■ 

IT) . 

When cosmic-ray protons and nuclei enter the atmosphere, they interact and pro- 
duce all kinds of secondary particles, which in turn interact or decay or propagate to 
the ground, depending on their intrinsic properties and energies. In the GeV energy 
range, the most abundant particles at the ground are the neutrinos. At produc- 
tion, there are approximately twice as many muon neutrinos as electron neutrinos. 
This characteristic ratio arises from the pion-muon-electron decay chain in which 
a muon and one muon-neutrino (or anti-neutrino) are produced when the charged 
pion decays, while the subsequent muon decay produces an electron neutrino as well 
5_i ■ as another muon neutrino. 

In 1960 Markov P suggested using Cherenkov light in a lake or the deep ocean 
to do neutrino physics with the atmospheric neutrino beam; in particular, to inves- 
tigate the question whether electron and muon neutrinos are distinct species. About 
the same time Greisen [2| described a proposal to search for astrophysical sources 
of neutrinos with a water Cherenkov detector in a deep mine. Atmospheric neutri- 
nos were, however, first detected with electronic detectors-as horizontal neutrino- 
induced muons in a mine in South Africa ,3J and in an iron calorimeter in a deep 
mine in the Kolar Gold Fields of India [4 . More than a decade passed before a large- 
scale effort to instrument a large body of water was undertaken by the DUMAND 
Project 

Atmospheric neutrino fluxes were first calculated in the 1960s. One approach jB] 
is to infer the neutrino fluxes from measurements of closely related muons. More 
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recently this approach has been refined to take account of effects of the geomagnetic 
field and of the separate contributions of pions and kaons [Zj . Most later calculations 
follow Refs. [HE] and calculate both the muon [8J and the neutrino fluxes [9] starting 
from the primary cosmic-ray spectrum at the top of the atmosphere. 

With the advent in the 1980s of large underground detectors to search for proton 
decay it became possible to measure increasingly large samples of events induced 
by neutrinos of both types (and from all directions because the Earth is transparent 
to neutrinos). Such events were of interest primarily as the background for proton 
decay. Both water Cherenkov detectors [TU] ITT] and segmented iron calorimeters j3J 
UU EU were used. 

Evidence gradually accumulated that the ratio of the two neutrino flavors was 
different from the expected value of two. The first hint came from the 1MB experi- 
ment [T3] , which reported too few muon decays compared to what was expected from 
from interactions of inside their detector. The 1MB experiment later concluded, 
however, that their measurement of the ratio of stopping to throughgoing muons 
was inconsistent with oscillations 16 . The Kamiokande experiment [Tj\ found an 
anomalously low ratio of / v e and suggested neutrino oscillations as a possible 
explanation. As the 1MB experiment accumulated more data JH] , they too found a 
persistently low j v e flavor ratio. Meanwhile, however, the Frejus experiment JH] 
measured a ratio consistent with the expectation. After some six years of running 
with the shielded Kamiokande II-III detector, enough data was accumulated to see 
a suggestion of the pathlength dependence expected from oscillations with parame- 
ters such that upward z/ M in the multi-GeV energy region oscillate while downward 
neutrinos do not [T9] . 

By measuring the directions and energies of thousands of neutrinos of both fla- 
vors in the past decade, the Super-Kamiokande Collaboration has shown that the 
anomalous ratio is a consequence of neutrino flavor oscillations during their propaga- 
tion from the atmosphere to the detector [20J. They have measured the oscillations 
in the muon neutrino-tau neutrino sector [21] by fitting the energy and angular 
dependence and other observed features of their data, while ruling out oscillation 
to sterile neutrinos at 99% confidence level. It is now possible to demonstrate the 
L/E pathlength/energy dependence expected for oscillations [22] ■ Best fit param- 
eters in a two-flavor approximation [y^ <-> v T ) are 5m 2 = 2.1 x 10~ 3 eV 2 and 
maximal mixing [2H]. Results from MACRO [21] and Soudan are consistent 
with these parameters. Measurements of atmospheric neutrinos are reviewed in 

Refs. [23 123123 • 

With the parallel discovery [2H1 I2H1 ISO] of oscillations of solar electron-neutrinos 
in a different region of parameter space, a pattern of mixing among the three neu- 
trino types (electron, muon and tau) is beginning to emerge in which two mixing 
angles are large while one is smaller (though not yet measured). Eventually the na- 
ture of the full neutrino mass matrix will be addressed with long-baseline accelerator 
neutrino beams by experiments now under construction or planned. This will take 
some time, however. Accordingly, there is interest in refining calculations of the 
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atmospheric neutrino beam to the point where effects of sub-dominant mixing may 
be resolved. To take full advantage of the power of the Super-Kamiokande experi- 
ment requires reducing uncertainties in the absolute normalization of the neutrino 
flux and in the ratios of electron to muon neutrinos and neutrinos to anti-neutrinos. 
The normalization depends on the primary cosmic-ray intensity and on the details of 
production of pions and kaons in cosmic-ray interactions in the atmosphere. The ra- 
tios depend primarily on pion and kaon production. This paper reviews the present 
level of uncertainties. 

2 Calculation 

The neutrino flux is a convolution of the primary cosmic-ray spectrum, modified by 
geomagnetic cutoffs, with the yield of neutrinos per incident cosmic-ray. Schemati- 
cally, 

<p u = <p p ® R p ® Y p ^ u + <Pa ® Ra <8> Ya-^v (1) 

A 

The first term represents the contribution of free protons and the second term the 
contribution of nucleons bound in nuclei. The two need to be calculated separately 
even in the superposition approximation in which it is assumed that bound nucleons 
interact independently as if they were free. This is because the geomagnetic cutoff 
depends on rigidity (total momentum per unit charge) whereas pion production 
depends on energy per nucleon. In each term the factor (j){E) represents the primary 
cosmic-ray spectrum as a function of energy per nucleon. R(E, 0®, 0®, 9, 0) indicates 
the cutoff rigidity (more exactly the transmission probability) which depends on 
latitude and longitude and on the local zenith and azimuth of the arrival direction 
of each primary cosmic ray. The yield Y(Ea, E u ) gives the number of neutrinos with 
energy E v produced per primary of energy Ea and is obtained by calculating the 
cosmic-ray induced cascades in the atmosphere. A complete review of the calculation 
of the flux of atmospheric neutrinos is given in Ref. In this paper I illustrate 
the considerations involved in calculating the production spectrum of atmospheric 
neutrinos by comparing results of three recent calculations (22 Elj . 

All these calculations are three-dimensional, taking account of the deviations 
of the neutrinos from the directions of the primaries that produce them. In ad- 
dition, the calculations of Refs. [33] and j^lj also account for the effect bending 
of muons in the geomagnetic field before they decay. There are several other fully 
three-dimensional calculations, including the recent papers of Refs. jHSHSHlEZj which 
appeared after the review [31]. The challenge of the full three-dimensional calcu- 
lation is indicated by the five arguments of the cutoff function above. Cascades 
must be generated for relevant primary energies all over the globe (#©,0®) taking 
account of the cutoff for primaries from all directions (0,0). Since most of the 
created neutrinos miss the detector, the calculation is highly inefficient. Although 
three-dimensional aspects of a calculation are technically challenging, they are not 
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the most important sources of uncertainty in the calculated fluxes of atmospheric 
neutrinos. The biggest uncertainties come from the primary cosmic-ray spectrum 
and the treatment of hadronic interactions. 



2.1 Primary spectrum 
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Figure 1: The flux of nucleons. The heavy black line shows the numerical form of 
Eq. El The lighter lines show extrapolations of fits [10] to measurements of protons, 
helium and three heavier groups below 100 GeV/nucleon. 



In the past decade new measurements of the primary spectrum with magnetic 
spectrometers have improved our knowledge of the primary spectrum. In particu- 
lar, the AMS [HB] and BESS [30] detectors give results for the proton spectra up 
to 100 GeV that agree with each other within 5%. As a consequence, current cal- 
culations of atmospheric neutrinos are using fits to the primary spectrum in which 
these data have been given priority. An example of such a fit [10] is shown by 
the thin lines in Fig. ^ Here the various nuclear groups are plotted as nucleons 
per GeV/nucleon, which is the quantity most directly related to secondary particle 
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production. Free protons make up about 70% of the flux of all nucleons, helium 
about 20%, and heavier nuclei the rest. Another recent series of spectrometer mea- 
surements (CAPRICE |41J) continues to give results for protons about 15% lower 
than BESS and AMS, which can be taken as an indication of the uncertainty in the 
absolute normalization of this component of the primary spectrum. It should also 
be noted that BESS [HE] and AMS [12] results for helium do not agree with each 
other as well as for protons. 

Measurements with spectrometers presently extend up to about 500 GeV only [33] . 
Higher energy data are from balloon-borne ionization calorimeters, which sample the 
fraction of energy deposited in electromagnetic cascades inside the calorimeter. Sys- 
tematic errors in assigning primary energy are larger as a consequence. Data from 
JACEE jS] and RUN JOB EE] are included in Fig. [TJ Preliminary data from the 
ATIC experiment |J7| fill in the gap in the 1-100 TeV region. The ATIC data for 
protons appear consistent with all previous data. Their helium data, if confirmed, 
indicate a preference for the lower RUNJOB normalization for helium. 
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Figure 2: Ratio of the all-nucleon fluxes used in several calculations to the power-law 
form of Eq. El 

The heavy solid line in Fig. ^ is a simple power-law approximation to the spec- 
trum of all nucleons, 
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where E is total energy per nucleon. Fig. [21 displays the all-nucleon spectra used to 
calculate the neutrino flux as a ratio to the power-law form of Eq. [21 Refs. [321 EH! 
use the fits of Ref. [HI] 5 while Ref. [33| use slightly different parameters for the 
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heavy components at low energy and a slightly flatter power-law extrapolation (- 
2.71) above 100 GeV. Nevertheless, the primary all-nucleon spectrum of the three 
calculations [321 ESI E3] are the same to within about 1%. Moreover, simple power 
law approximation |21 to the all-nucleon spectrum represents these fits to better than 
10% over the whole energy region. For reference, Fig. |2] also shows the primary 
spectrum used in an earlier calculation [1H], which is somewhat lower than the new 
fits below 100 GeV and significantly higher at high energy. 



Comparison of three neutrino fiux calculations at Super-K 
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Figure 3: The production spectrum of + and v e + v e integrated over the 
atmosphere and averaged over all directions at Super-K (left) and Soudan or SNO 
(right). 



2.2 Treatment of hadronic interactions 

The representations of hadronic interactions used in Refs. |S2l EI EH] are completely 
independent. As a consequence, it is not surprising that they give different results. 
Because the assumptions about primary spectrum are essentially identical, differ- 
ences among the results of the three calculations reflect the level of uncertainty due 
to treatment of hadronic interactions. For low energies the fluxes depend also on the 
geomagnetic cutoffs, so comparison must be made for the same location. Figure El 
shows the atmospheric neutrino fluxes from the three calculations [S2l ESI El] • Dif- 
ferences among the calculations are at the level of 10% at Super-K. The differences 
are somewhat larger for the sites at high geomagnetic latitude (where the cutoff 
for downward primaries is negligible) , which probably indicates that the low-energy 
hadronic interactions are less well understood. 

The flavor ratios for the three calculations are shown for the low-energy region 
in Fig. |3J along with the neutrino/antineutrino ratios. Differences are at the level 
of 3% in this energy range for the flavor ratios and 5% for the v / v ratios. 
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Figure 4: The ratio + to v e + u e (left panel) and vjv (right panel) at Super-K. 
The production spectra are integrated over the atmosphere neglecting oscillations 
and averaged over all directions before taking the ratio. 
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Figure 5: Distribution of neutrino energies (left) and primary energy per nucleon 
(right) giving rise to various classes of events. 



3 Higher energies 

Response functions for various measurements of atmospheric neutrinos are given in 
Ref. [21] and reproduced here in Fig. Sub-GeV and multi-GeV events involve 
neutrinos with energies below 2 GeV and from 1 to 20 GeV respectively. Upward 
stopping muons are produced by muon neutrinos with energies in the range 3-100 
GeV, while upward throughgoing muons come from neutrinos with energies in the 
range 10 GeV to 10 TeV. The relevant primary energies per nucleon are roughly a 
factor of ten higher, but with rather broader distributions. Thus the normalization of 
the sub-GeV neutrinos is determined almost entirely by primaries with energies less 
than 100 GeV/nucleon and hence covered by the most precise magnetic spectrometer 
measurements. On the other hand, the primary response function for neutrino- 
induced muons extends to 10 TeV and beyond and hence is subject to somewhat 
larger uncertainty in normalization. The difference in response function for neutrino- 
induced upward muons at Super-K from that at AMANDA arises from the higher 
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muon energy threshold in AMANDA, which is of order 100 GeV at the interaction 
vertex [4*^] . 



Comparison of v u +v u fluxes 
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Figure 6: Comparison of fluxes of muon neutrinos at high energy. 

Again we can make use of the fact that the assumed primary spectra in the 
three calculations EH1 EH] are almost the same to illustrate differences due to the 
treatment of hadronic interactions. Fig. |U] compares fluxes of muon neutrinos for 
these three calculations in the energy range responsible for neutrino-induced muons. 



3.1 Analytic approximations 



At high energy, numerical integration of the analytic cascade equations gives suf- 
ficiently accurate results to be useful for understanding the characteristic features 
of the results. The asterisks in Fig. |H1 show the numerical results of Lipari |3U| , 
including the contribution from decay of muons. 1 At high energy most muons reach 
the ground before decaying and therefore do not contribute to the neutrino flux. 
For E v > 100 GeV for example, less than 15% of muon neutrinos are from decay of 
muons, and the fraction decreases further as energy increases. 

Neglecting neutrinos from muon decay, a simple approximation for the flux of 



+ Ufj, from decay of pions and kaons is 
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lr rhe primary spectrum of |5l)j has the same form as Eq. [3 with a higher normalization. 
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where 4>n(E u ) is the primary spectrum of nucleons evaluated at the energy of the 
neutrino (see Eq. |5J). The constants = m^/m| for i = (tt,K), while the constants 
B depend on the hadron attenuation lengths as well as decay kinematics. The 
critical energy for pions is e n ~ 115 GeV, while for kaons e# ~ 850 GeV [oTj . 
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Figure 7: Spectrum-weighted Z-factor for K + production by nucleons. 

Asymptotically, the neutrino spectrum is one power steeper than the primary 
spectrum. For fixed zenith angle 9 the transition occurs first for neutrinos from pion 
decay because of its lower critical energy. As a consequence, the relative contribution 
from kaons increases with energy. Kaons become the dominant source of neutrinos 
for E v > 100 GeV. (The onset of kaon dominance is somewhat slower at large angle 
because of the cosine factor in the denominator of the kaon term in Eq.EJ) Important 
differences among the calculations can be traced to differences in treatment of kaon 
production, in particular to the channel p — > AK + . Fig. [7| shows the spectrum- 
weighted moments for four calculations. 

Several features of Fig. |H1 can be understood with the help of Eq. El 

• The general decrease of the vertical to horizontal ratio at large energy, which 
appears in all calculations, is a consequence of the factor of cos(#) in the 
denominators. 
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• The primary spectrum of [31] is significantly lower than that of [48 at high 
energy, as shown in Fig. [21 The calculation of Ref. [31] is nevertheless high, 
comparable to that of Ref. [2H] because the kaon production assumed in the 
new calculation is larger, as shown in Fig. [7| 

• The ratio (> 1) of the flux of Ref. [31] compared to that of Ref. [33]; and its 
increase with energy, may also be attributed to the differences in the assump- 
tions for kaon production as shown in Fig. [7| 

3.2 Muon fluxes 

Compilation of muon spectra 
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Figure 8: Summary of measurements of the vertical muon intensity at the ground. 
The solid line shows an analytic calculation [52 . The dotted line shows the spec- 
trum in the absence of decay and energy loss, or equivalently the muon production 
spectrum integrated over the atmosphere. 

Calculation of the flux of atmospheric muons follows the same form as Eq. Q 
with the yield of neutrinos replaced by the yield of muons. Because of the close 
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relation between muons and neutrinos, comparison of its predicted muon intensity 
to measurements is an important check on any calculation of atmospheric muons. 
Comparison to muon fluxes also serves to check the normalization of the calcula- 
tion. Because muons lose about 2 GeV in passing through the atmosphere, a ground 
level measurement probes high energies and is therefore most relevant to multi-GeV 
neutrinos and to neutrino-induced upward muons. Because of the kinematics of the 
7i — > fi decay, however, pion decay remains the dominant source of muons at all en- 
ergies. As a consequence, the leverage of the muon measurements for controlling the 
neutrino flux is limited above ~ 100 GeV where neutrino production is dominated 
by the contribution from kaons. 

Fig. |H1 contains a summary of measurements of the spectrum of vertical muons at 
the ground. All measurements have been corrected to sea level. At very high energy 
the muon spectrum becomes one power steeper than the parent spectrum of nucleons 
as a consequence of the extra power of 1 / E n in the ratio of pion decay length to 
interaction length, which reflects the decreasing probability of decay relative to re- 
interaction for charged pions at high energy. For E < e n ~ 115 GeV essentially 
all pions decay, and the muon production spectrum has the same power behavior 
as the parent pion and grandparent nucleon spectrum (a » 2.7). At low energy, 
however, muon energy- loss and decay become important, and the muon spectrum at 
the ground falls increasingly below the production spectrum. To account for all the 
complications one generally resorts to Monte Carlo calculations. However, analytic 
approximations of the effects are also possible |50| . Fig. |H1 shows a comparison of one 
such calculation which uses as input the simple power- law primary spectrum 
of Eq. [21 Data in Fig. |H] are from many measurements with spectrometers at the 
surface [EH EH ESI EH E2 EE] and underground [SU ED] . 

The comparison of the analytic calculation with the data in Fig. |S] leads to two 
remarks. First the good agreement gives a general confirmation of the primary 
spectrum at high energy. Second, it serves to illustrate the level of systematic 
differences among the various measurements. 

4 Concluding remarks 

In their recent detailed paper, [2H] the Super-Kamiokande group give an extensive 
discussion of how the best fit oscillation parameters are determined from a com- 
parison of their data with simulations based on the atmospheric neutrino flux of 
Ref. (331- In particular, their Table VII enumerates the assumed uncertainties in 
various properties of the calculated neutrino flux and the shifts in those properties 
needed to obtain the best fit. Apart from the overall increase in normalization, 
the largest shifts relative to what is assumed in Ref. [33] concern the slope of the 
primary spectrum at high energy and the K/pi ratio. The latter is decreased by 
6% overall, while the primary spectrum is significantly harder as compared to the 
fits of Ref. |1D] (from -2.71 to -2.66 above 100 GeV). As noted in the discussion of 
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Figs. El and the high energy t/-spectrum could also be increased by enhancing kaon 
production preferentially at high energy. 

The Super-K fitting procedure uses the atmospheric neutrino data itself to cor- 
rect the calculated neutrino flux. The neutrino data have the advantage in principle 
that (apart from oscillations!) the intensity observed at the detector integrates over 
the whole atmosphere without complications of energy loss and decay that dom- 
inate the observed muon flux. For example the required increase in the overall 
normalization of 11.9% could be interpreted as favoring the higher normalization of 
Refs. (2HI GUI] as compared to that of Ref. [JT]. The observations depend, however, 
on properties of neutrino interactions as well as on the neutrino flux. As properties 
of neutrino interactions become better determined, for example from near detectors 
of long-baseline experiments, this iterative procedure may become an important way 
to refine our understanding of the atmospheric neutrino flux at production. 

Acknowledgments. I am grateful to G. Battistoni and M. Honda for providing 
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